High-Order Coupled Cluster Calculations Via Parallel Processing: 

An Illustration For CaV 4 9 
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The coupled cluster method (CCM) is a method of quantum many-body theory that may provide 
accurate results for the ground-state properties of lattice quantum spin systems even in the presence 
of strong frustration and for lattices of arbitrary spatial dimensionality. Here we present a significant 
extension of the method by introducing a new approach that allows an efficient parallelization of 
computer codes that carry out "high-order" CCM calculations. We find that we are able to extend 
such CCM calculations by an order of magnitude higher than ever before utilized in a high-order 
CCM calculation for an antiferromagnet. Furthermore, we use only a relatively modest number 
of processors, namely, eight. Such very high-order CCM calculations are possible only by using 
such a parallelized approach. An illustration of the new approach is presented for the ground-state 
properties of a highly frustrated two-dimensional magnetic material, CaV-iOg. Our best results for 
the ground-state energy and sublattice magnetization for the pure nearest-neighbor model are given 
by Eg/N — —0.5534 and M = 0.19, respectively, and we predict that there is no Neel ordering in 
the region 0.2 < J2/J1 < 0.7. These results are shown to be in excellent agreement with the best 
results of other approximate methods. 
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A new procedure for carrying out high-order coupled 
cluster method (CCM) [1IIB&I1I1lEI3 calculations 
via parallel processing is presented here. The CCM may 
be applied to systems demonstrating strong frustration 
and for arbitrary spatial dimension of the lattice. We 
illustrate our new app roach by applying it to CaV/iOg 
HEEimHIlim at zer0 temperature. An in- 
creasing number of insulating quantum magnetic systems 
for lattices of low spatial dimensionality are being stud- 
ied experimentally. Indeed, the calcium vanadium oxide 
(CAVO) materials are one particularly useful example. 
They exhibit strong frustration for a number of different 
crystallographic lattices with respect to varying chemical 
composition, and they may demonstrate "novel" ground- 
state ordering. Indeed, theoretical evidence suggests that 
CAVO systems may contain a number of differing quan- 
tum ground states (see, e.g., Refs. HEI) at zer0 tem- 
perature as a function of varying nearest-neighbor next- 
nearest-neighbor bond strengths. The relevant Hamilto- 
nian is given by 
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where i runs over all lattice sites, and j and k run over 
all nearest-neighbor and next-nearest-neighbor sites to 
i, respectively, counting each bond once and once only. 
The bond strengths are given by J\ and J2 for nearest- 
neighbor and next-nearest-neighbor terms, respectively. 
The lattice and exchange "bonds" are illustrated graph- 
ically in Fig. 

For this model, collinear Neel ordering is observed with 
respect to nearest neighbors at J\ > and J2 = and 
with respect to next-nearest neighbors at J2 > and 
J\ = 0. Exact diagonalizations predict an interme- 
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FIG. 1: The unit cell of this system contains four sites and is 
indicated in the figure. Nearest-neighbor bonds are shown by 
the full lines and next-nearest-neighbor bonds by the dashed 
lines. 



diate phase (or phases) for 0.2 < J2/J1 < 0.7, which 
is also in good agreement with results of a "spin-wave 
theory-like" treatment that predict such an inter- 
mediate regime for 0.25 < J%j J\ < 0.8. Strong evi- 
dence has been gathered using approximate approaches 
H El- that the ground state in this "interme- 

diate" regime demonstrate dimer or plaquette ordering. 
However, conclusive evidence for the nature of the ground 
state across the whole of this regime remains to be yet de- 
termined. The ground state of the pure nearest-neighbor 
model at J = 1 and J2 = demonstrates semi-classical 
Neel order with a sublattice magnetization of 0.178(8) 
via QMC [13 ■ 

We now turn our attention to the application of the 
CCM to this model. We choose two model states that 
are the classical ground states of this model in the lim- 
its J2/J1 — » and J2/J1 — * 00 (with both J2 and 
Ji positive). These are, namely, one in which nearest- 
neighbor spins are antiparallel, henceforth denoted the 
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FIG. 2: CCM results for the ground-state energy per spin at 
Ji — 1 and J2 — using the LSUBm approximation with 
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FIG. 3: CCM results for the ground-state energy per spin of 
the J1-J2 model (with Ji = 1) using the LSUBm approxi- 
mation with m = {6, 8, 10}. We see that our results converge 
rapidly with m across a wide range of J2 up to and including 
the point J2 = 0. Extrapolations using Pade approximants 
are also shown. 



nearest-neighbor Neel model state; and a second model 
state in which next-nearest-neighbor spins are antipar- 
allel, henceforth denoted the next-nearest-neighbor Neel 
model state. We use four spins in the crystallographic 
primitive cell in both cases. This unit cell is indicated 
in Fig. We perform a notational rotation of the "up" 
spins to "down" spins, as in Refs. @, II 11 H H II 11 
for example. We use the full symmetries of the lattice 
to further reduce the computational problem of finding 
the fundamental CCM configurations and in determining 
and solving the CCM equations in both cases. 

The exact ket and bra ground-state energy eigenvec- 
tors, I*) and (\E'|, of a general many-body system de- 
scribed by a Hamiltonian H , 

H\*) = E g \V); (*\H = E g {*\, (2) 

are parameterized within the single-reference CCM as fol- 



lows: 

|*)=e s |$) ; S = Y,SiC+, 

<*IH$|£e- s ; s = i + J2SiC7- (3) 

The single model or reference state |$) is required to have 
the property of being a cyclic vector with respect to two 
well-defined Abelian subalgebras of multi-configurational 
creation operators {C^} and their Hermitian-adjoint de- 
struction counterparts {CJ = (Cj)'}. The determina- 
tion of the correlation coefficients {Si, Si} is achieved by 
by requiring the ground-state energy expectation func- 
tional H({Si, Si}), to be stationary with respect to vari- 
ations in each of the (independent) variables of the full 
set. We thereby easily derive the following coupled set of 
equations, 

SH/5Si = => (<I>\C]-e- s He s \<5>} = 0, 1^0 ; (4) 
SH/SSi = => (<S>\Se- s [H, C+]e s |$) =0, I ^ (5) 

Equation (0J also shows that the ground-state energy at 
the stationary point has the form 

E g = E g {{Si}) = (*|e- s fre 5 |*> . (6) 

It is important to realize that this (bi-)variational for- 
mulation does not lead to an upper bound for E g when 
the summations for S and S in Eq. J3J are truncated, 
due to the lack of exact Hermiticity when such approx- 
imations are made. However, one can prove that the 
important Hellmann-Feynman theorem is preserved in 
all such approximations. We utilize various approxima- 
tion schemes within S and S. The three most commonly 
employed schemes previously utilized have been: (1) the 
SUBn scheme, in which all correlations involving only 
n or fewer spins are retained, but no further restriction 
is made concerning their spatial separation on the lat- 
tice; (2) the SUBn-m sub-approximation, in which all 
SUBn correlations spanning a range of no more than m 
adjacent lattice sites are retained; and (3) the localized 
LSUBm scheme, in which all multi-spin correlations over 
all distinct locales on the lattice defined by m or fewer 
contiguous sites are retained. 

We now solve the CCM equations in parallel by re- 
arranging Eqs. and © for the ket and bra states, 
respectively, where 

Si = fi( Si, Si-i, Si + i, Sftf) ; (7) 
Si = fi( Si, ■ ■ ■ , Sfff) 

Si, S1-1, Si + i, SN f ) ; (8) 

and where Nf is the number of fundamental configu- 
rations. There are therefore Nf equations for both // 
and //, which refer to the ket and bra states, respec- 
tively. Each equation contains a finite-number of terms 
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FIG. 4: CCM results for the sublattice magnetization of the 
J1-J2 model using the LSUBm approximation with m = 
{6,8,10}. Extrapolations using Pade approximants are also 
shown. 



TABLE I: CCM LSUBm results for the CAVO antiferromag- 
net. The numbers of fundamental configurations nearest- 
neighbor and next-nearest-neighbor Neel model states are 
given by Nf' n ' and AfJ'"'"', respectively. Results for the 
points at which the sublattice magnetization goes to zero for 
the nearest-neighbor and next-nearest-neighbor Neel model 
states are indicated by J2/J1I1 and J2/J1I2, respectively. 
(Errors for the critical points are in the last decimal place 
shown.) Ground-state energies and sublattice magnetiza- 
tions for nearest-neighbor bonds only, namely at Ji = 1 
and J2 = 0, are extrapolated using Pade approximants for 
m — {4, 6, 8, 10}, where I is the order of the polynomial in the 
denominator. By way of comparison, results of exact diago- 
nalizations (ED) for E g /N and M of Ref. 0] and values for 
J2/J1I1 and J2/J1I2 of Ref. 0] are also given. 
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and that, indeed, this is always the case for the CCM 
for Hamiltonians that contain a finite number of creation 
and destruction operators. Our approach is to iterate 
directly the set of equations in (JJJ to convergence first 
in order to obtain the ket-state correlation coefficients. 
We then iterate directly the set of equations in JSJ) to 
convergence in order to obtain values for the bra-state 
correlation coefficients. 

Our parallelization technique is now to split the prob- 
lem of determining and solving the Nf equations for both 



the bra and ket states of Eqs. Q and (|SJ) between each 
processor in the parallel cluster. Thus, each equation for 
I = {1, • • • , Nf} is ever dealt with by one node in the 
cluster only and the results of each node are collected 
together at each iteration. (Again, we remember that 
we iterate the ket-state equations first and then the bra- 
state equations.) Thus, the set of equations for the ket 
and bra states are solved independently and in parallel. 
Eight Compaq "Wildfire" processors were used here, each 
with 2 GIGs of RAM and CPU speeds of 731 MHz. The 
run for LSUB10 with the nearest-neighbor model state 
with 22395 fundamental configurations took 11 days for 
18 distinct values of Jij J\. 

We now wish only to illustrate this approach by consid- 
ering the the CAVO system of Eq. QJ. The ground-state 
energies of the pure nearest-neighbor antiferromagnet (at 
J\ = 1 and J2 = 0.0) are shown in Table [I] Fig. [3 shows 
LSUBm ground-state energy plotted against rrT 2 for 
rn = {4, 6, 8, 10}. "Linear extrapolations" of data points 
have previously been used to great effect for isotropic 
nearest- neighbor antiferromagnets 0, EJ 0, El • In this 
case of the CAVO model studies here, we obtain a value 
of Eg/N = —0.55310(5). We also carry out extrapo- 
lations using Pade approximants 01 for E g /N plotted 
against mT 2 for m = {4,6,8,10}. I is the order of the 
polynomial in the denominator and that we perform an 
exact fit of the curve to the data. Hence, in contrast 
to the simple linear fit, there is no "curve-fitting" er- 
ror. Extrapolations using Pade approximants for the 
ground-state energy at J\ = 1 and J2 = 0.0 are found 
to be in excellent agreement with a variational QMC re- 
sult of E g /N = -0.5501 and an ED result [3 of 
Eg/N = —0.5533. We conclude that our results are con- 
verged to at least 3 decimal places, and we believe that 
these results constitute one of the most accurate esti- 
mations of the ground-state energy for this system ever 
found. 

The number of fundamental configurations for the 
SUB10-10 approximations were 22340 and 39434 for 
the nearest-neighbor and next-nearest-neighbor model 
states, respectively, and this calculation for the next- 
nearest-neighbor model state is the largest CCM calcula- 
tion as yet performed for quantum antiferromagnetic spin 
systems by an order of magnitude compared to the previ- 
ous highest . The results for the ground-state energies 
for both model states are shown in Fig. |2]with respect to 
J2/J1 ■ We see that these results converge quickly with 
truncation index m for the LSUBm approximation over 
a wide range of Ji/ J\- The wide area without results is 
an indication ground-state ordering (e.g., dimer or pla- 
quette) not accessible with our Neel model states. 

LSUBm results for the sublattice magnetization at 
J\ = 1 and J2 = are again shown in Table [I] In 
contrast to the case for the ground-state energy, extrap- 
olations for the sublattice magnetization do not seem to 
follow a similar linear extrapolation "rule" as shown in 
Fig. 0] However, as in Refs. @, 0, 0, IE Hi we may plot 
the sublattice magnetization against m _1 and use Pade 
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approximants in order to carry out the extrapolations. 
This gives values for M in the range 0.19 to 0.20, which 
is again in good agreement with the best of other approx- 
imate methods, such as results from QMC [lOj and ED 
E9 of M = 0.178(8) and M = 0.2303, respectively. We 
therefore use Pade approximants for increasing values of 
J2/J1 > as a reliable indicator of the position of the 
phase transition points is given by the position at which 
the sublattice magnetization goes to zero as a function of 
J2/J1. Fig. 0| plots our high-order CCM results for the 
amount of sublattice magnetization. We may see from 
Table [Q that our results for the positions of the phase 
transition points are in good agreement with results of ex- 
act diagonalizations [uj that predict that plaquette and 
dimer ordering exists in the region 0.2 < J2/J1 < 0.7. 

We have shown in this article that it is possible to 
simulate the properties of a CAVO material at zero tem- 
perature using parallelized high-order CCM techniques. 
The construction of CCM codes to such high orders of 
approximation is a non-trivial task. Furthermore, the 
parallelization of the CCM code adds an additional level 
of complexity to it. The 39434 fundamental configura- 
tions for the next-nearest-neighbor Neel model state is 
a massive leap forward in the practical application of 
the CCM, which is impossible without an efficient paral- 
lelization. However, the numbers of fundamental config- 
urations of order 10 6 for quantum antifcrromagnets are 
clearly within our reach via parallel processing for even 



larger clusters of processors than the relatively modest 
number of (eight) processors used here. We have there- 
fore "proven the principle" that the CCM can be ef- 
fectively parallelized, and we believe that we have also 
shown here that this, in turn, leads to great increases in 
accuracy of observed expectation values. 

The extension of these calculations to treat quantum 
systems that demonstrate "novel ordering," such as with 
plaquette or dimer solid ground states, is, in principle, 
straightforward. The basic building blocks for the CCM 
using dimer models states are already in place, and low- 
order calculations have already been attempted success- 
fully for a one-dimensional frustrated spin system 
However, procedurally, the spatial dimension of the lat- 
tice makes no difference to the coding of the problem, 
although the computational expense grows with it and 
so efficient parallelized approaches, such as that shown 
here, would be very useful in this case. Extensions of the 
method to consider excitation energies have been very 
successful 0, and it is again a straightforward matter 
to obtain excitation spectra using "localized" high-order 
CCM approximations. The CCM may be applied to con- 
sider systems even in the presence of strong frustration 
and for lattices of arbitrary spatial dimensionality, and 
the CCM now constitutes a very powerful and a useful 
tool in order to understand the basic properties of quan- 
tum spin systems. 
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